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ABSTRACT 


Technical  Report 

This  report  presents  roll  damping  coefficients  and  pitch  damping  coefficients 
obtained  from  dynamic  rolling  tests  and  static  wind  tunnel  tests  of  the 
Australian  81mm  Improved  Mortar  Projectile,  IMP.  An  80%  scale  model  was 
used  in  the  dynamic  roil  tests  and  a  full  scale  model  was  used  in  the  static  wind 
tunnel  tests. 
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Nomenclature 


Ci  rolling  moment  coefficient,  torque  coefficient  I 

C|0  static  rolling  moment  coefficient 

Qp  roll  damping  coefficient 

C,,,^  pitch  damping  coefficient 

D  drag  force  I 

1  rolling  moment 

L  lift  force 

m  pitching  moment 

X  axial  force  ^ 

Z  normal  force 

0  pitch  angle  (deg) 
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1.  Introduction 

The  Australian  Army's  Engineering  Development  Establishment  (EDE)  is  working  with 
Australian  Defence  Industries  Pty  Ltd  (ADI)  to  develop  and  produce  a  new  improved  81mm 
mortar  projectile  to  replace  the  M374  High  Explosive  (HE)  round.  The  Weapon  Aerodynamics 
Discipline  of  the  Aeronautical  Research  Laboratory  (ARL)  was  requested  by  Explosive 
Ordnance  Division  (EOD)  of  Materials  Research  Laboratory  (MRL)  to  determine  pitch  and 
roll  damping  data  (Q^  and  C\p  respectively)  for  the  81mm  mortar  from  wind  tunnel  tests.  The 
data  will  be  added  to  existing  static  aerodynamic  data  obtained  from  previous  wind  tunnel  tests 
of  an  80%  and  100%  scale  mortar  model  (references  1  and  2),  and  will  be  used  to  model  the 
flight  behaviour  of  the  mortar  and  to  establish  its  stability  and  range  and  the  effect  of  launch 
transients. 


2.  Test  Details 

2.1  Wind  Tunnel 

The  wind  tunnel  used  for  these  tests  was  the  ARL-Salisbury  SI  wind  tunnel  which  is  a  closed 
circuit  continuous  operation  tunnel  The  working  section  has  dimensions  of  380  mm  x  360  mm 
with  slots  and  has  the  capability  for  Mach  numbers  of  0.35  to  1.0  and  1.4  to  2.8.  Tunnel  flow 
conditions  are  set  and  recorded  using  a  static  pressure  port  upstream  of  the  model,  a  pitot  tube 
located  upstream  of  the  contraction,  and  a  temperature  probe  in  the  settling  chamber.  The 
model  was  positioned  in  the  working  section  by  a  pitch  and  roll  mechanism  located  beneath  the 
working  section. 

2.2  Models 

Models  tested  were  an  80%  scale  metal  model  and  a  full  scale  metal  model  Doth  models  were 
designed  and  manufactured  at  ARL,  based  on  drawings  of  the  full  scale  Improved  Mortar 
Projectile  (IMP)  supplied  by  EDE  and  ADI  80%  scale  and  full  scale  HL  18622-022  fuzes  and 
DE  132410018  (extruded,  canted  metal)  fins  were  attached  to  the  corresponding  scale  model. 


2.3  Experimental  Technique 
2.3.1  Roll  Damping 

Roll  damping  coefficients,  C^,  were  obtained  from  dynamic  wind  tunnel  tests  carried  out  with 
an  80%  mortar  model  attached  to  a  roll  balance  as  shown  in  figures  1  and  2. 

The  roli  balance  is  fitted  with  a  vnicli  consists  of  an  expanding  fibrous  ring  operated  by 

compressed  air  working  on  pistons  in  the  balance  shaft  With  the  brake  on,  the  wind  tunnel 
was  started,  and  once  the  desired  Mach  number  was  reached,  the  brake  was  released,  allowing 
the  model  to  accelerate  from  rest  to  a  constant  roll  rate  During  this  lime,  roll  rates  and 
corresponding  times  were  recorded  and  saved  to  disk.  Tests  were  carried  out  at  Mach  numbers 
of  0.5,  0.7,  0.8,  0.9  and  0.95,  which  encompasses  the  flight  speed  range  of  the  81mm  IMP 
mortar,  at  pitch  angles  of0°,  5°  and  10°. 

Values  of  C|p  were  calculated  from  the  dynamic  rolling  data .  combined  with  static  rolling 
moment  daia  obtained  from  reference  1  A  detailed  descrintum  of  the  derivaaon  is  given  in 
Appendix  A. 
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2.3.2  Pitch  Damping 

Pitch  damping  coefficients  were  calculated  using  the  technique  given  in  reference  4.  This  uses 
static  aerodynamic  lift  data  for  lifting  surfaces  and  the  body  acting  separately  to  determine  the 
contribution  ol  each  to  the  overall  pitch  damping. 

The  fin  lift  characteristics  were  determined  from  wind  tunnel  tests  on  a  full  scale  model.  Tests 

were  conducted  at  Mach  numbers  of  0.5,  0.7,  0.8,  0.9  and  0.95  and  through  a  pitch  angle  range  » 

of  -5°  to  +5°,  first  on  the  model  with  fins  attached,  and  then  with  a  plain  cylindrical  tail  piece 

(no  fins)  attached.  The  difference  between  the  "fins  on"  and  the  "fins  off"  lift  curve  slopes  was 

taken  to  be  the  lift  curve  slope  of  the  fins  alone  operating  in  the  body  wake.  These  tests  were 

conducted  as  part  of  a  larger  test  program,  and  are  described  in  more  detail  in  reference  2. 

Appendix  C  describes  the  theoretically  based  derivation  of  the  body  lift  contribution,  and  the  ^ 

use  of  the  fin  and  body  contributions  to  calculate  the  overall  pitch  damping  coefficients. 


3.  Discussion  of  Results  » 

3.1  Roll  Damping 

Values  of  Cip,  derived  from  the  method  described  in  Appendix  A,  are  presented  in  table  1  and 
plotted  in  figure  3. 

> 


• 

Estimates  of  the  static  rolling  moment,  C|o,  can  be  obtained  from  the  dynamic  rolling  moment 
data  by  extrapolating  a  least  squares  straight  line  fit  to  die  data  back  to  zero  roll  rate  (sec 
figure  A3).  The  values  of  Qo  so  obtained  differ  from  C|o’s  reported  in  reference  I  by  up  to 
10%.  This  gives  a  useful  estimate  of  the  uncertainty  in  using  these  data  to  calculate  roll 

damping,  therefore  C)  values  stated  above  should  be  regarding  as  having  uncertainties  of  the  0 

order  of  ±10%. 

3.2  Pitch  Damping 

Values  of  derived  from  the  method  described  in  Appendix  C  are  presented  in  table  2  and 

plotted  m  figure  4  • 


Mach 

No. 

Roll  Damping  Coefficient,  C|p 

a  =  0°  a  -  5°  a  =  10° 

0.50 

-0.245 

-0.242 

-0.256 

0.70 

-0.259 

-0.256 

-0.281 

O.oO 

-0.268 

-0.260 

-0.302 

0.90 

-0.268 

-0.270 

-0.299 

0.95 

-0.300 

-0.3 14 

-0.347 

Table  1:  Roll  damping  coefficient,  Q 


Mach 

No. 

Cm, 

0.50 

-32.4 

0.70 

-39.7 

0.80 

-41.4 

0.90 

-49.3 

0.95 

-49.4 

Table  2:  Pitch  damping  coefficient, 


Figure  4  also  shows  values  of  for  an  M374  81mm  mortar  projectile  taken  from 
reference  3.  Although  this  is  a  different  projectile  shape,  the  similarities  are  sufficient  to  expect 
the  values  of  Cm  to  be  similar.  Comparison  of  the  two  sets  of  data  shows  that  for  Mach 
numbers  up  to  0.8,  reasonable  similarity  is  observed,  but  differences  of  practical  importance 
exist  at  the  higher  Mach  numbers.  The  divergence  of  results  above  Mach  0.8  may  be  due  to 
the  unusually  large  model  size,  resulting  in  significant  wind  tunnel  wall  interference  effects  at 
the  higher  Mach  numbers.  Reference  2  discusses  this  in  more  detail.  Because  of  these  effects, 
and  the  uncertainties  inherent  in  the  application  of  static  data  to  the  estimation  or  dynamic 
parameters,  the  pitch  damping  results  presented  here  are  estimated  to  have  uncertainties  of  up 
to  ±20%.  However,  for  a  fin-stabilised  ballistic  projectile  of  this  type,  this  level  of  uncertainty 
is  acceptable  in  terms  of  predicting  the  overall  flight  trajectory. 
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4.  Conclusion 

Roll  damping  coefficients  obtained  from  dynamic  wind  tunnel  tests  at  various  Mach  numbers 
and  pitch  angles  are  presented  in  table  1  and  figure  3.  These  data  are  estimated  to  contain 
uncertainties  of  ±10%  due  to  a  lack  of  low  roll  rate  data,  and  uncertainties  in  the  measurement 
of  roll  bearing  friction. 

Values  for  pitch  damping  coefficient  obtained  using  the  theory  presented  in  Appendix  V  (as 
obtained  from  reference  4)  and  static  wind  tunnel  results,  arc  presented  in  table  2  and  figure  4 
These  data  arc  estimated  to  contain  uncertainties  of  up  to  +20%  due  to  wind  tunnel  wall 
interference  and  uncertainties  inherent  in  the  method.  However,  from  comparison  with  other 

produces  acceptable 
conditions  considered 

here 


experiments,  it  appears  that  the  theory  presented  in  Appendix  C 
approximations  for  pitch  damping  coefficients  for  configurations  and 
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Figure  3:  Roll  damping  coefficient,  Cip,  vs  Mach  No. 
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Figure  4:  Pilch  damping  coefficient.  t’m  ,  vs  Mach  No. 


APPENDIX  A 


DERIVATION  OF  C|p 


Derivation  of  C|p 
Symbols  used: 

Ci  rolling  moment  coefficient 

Cj^  static  rolling  moment  coefficient 

C|p  roll  damping  coefficient 

d  reference  length  (diameter) 

p  roll  rate 

V  velocity 

The  data  for  .oil  rate  versus  time  (see  figure  A1  for  an  example  of  the  roil  rates  obtained),  were 
differentiated  to  produce  values  for  roll  acceleration  versus  time.  From  a  knowledge  of  the  roll 
inertia  of  the  model,  roll  accelerations  were  then  converted  to  roll  torques,  to  give  a  set  of  roll 
torques  versus  roll  rate.  Bearing  friction  torque,  as  calculated  in  Appendix  B,  was  then 
subtracted  from  the  roll  torque,  producing  values  for  aerodynamic  torque  versus  roll  rate.  The 
aerodynamic  torque  was  non-dimensionalised  to  obtain  rolling  moment  coefficients,  Q,  versus 
roll  rate. 


Roll  damping  coefficient,  Qp  was  calculated  using  the  following  equation: 


C'p 


'  8p  Id;' 


Due  to  the  model's  rapid  acceleration,  and  the  timcbase  on  the  Hewlett  Packard  digital 
oscilloscope  used  to  record  the  data  having  to  be  set  so  as  to  cater  for  the  high  roll  rates 
expected,  roll  rates  less  than  150  rad/s  (-  24  Hz)  could  not  be  recorded.  Due  to  the  lack  of  low 
roll  rate  data,  it  was  decided  to  include  zero  roll  rate  data  (static  rolling  moment  coefficients) 
measured  in  earlier  tests  and  reported  in  reference  1 .  These  data  did  not  agree  well  with  simple 
extrapolation  of  the  dynamic  data  (see  figure  A3),  throwing  some  doubt  on  the  intermediate  roll 
rate  data.  To  get  the  best  estimate  of  it  was  decided  to  use  only  the  two  most  reliable  Q 
values,  ie.  Ci  from  static  test  data  and  Ci  at  the  maximum  (constant)  roll  rate  Figure  A2 
shows  the  straight  line  drawn  for  Mach  0.95  using  only  these  two  extreme  va'ues.  The  gradient 
of  this  line  is  the  roll  damping  coefficient,  C|  ,  for  Mach  0.95.  C|p  values  for  all  other  Mach 
numbers  were  obtained  from  the  same  procedure. 

Note  that,  in  figure  A2  the  magnitude  of  the  bearing  friction  relative  to  the  aerodynamic  torque 
is  indicated  by  the  height  above  the  roll  rate  axis  at  which  maximum  roll  rate  occurs.  At  this 
steady  state  roll  rate,  bearing  friction  equals  aerodynamic  torque.  It  can  be  seen  that  maximum 
bearing  friction  is  about  5%  of  maximum  aerodynamic  torque,  and  hence  moderate 
uncertainties  in  the  measurement  of  bearing  friction  are  not  a  significant  contributor  to  the 
overall  uncertainty  of  the  result 
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Figure  A 1 :  Roll  rate  vs  time  (for  pitch  angle  -  0°) 
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Figure  A2:  Rolling  moment  coefficient  vs  roll  rate  (for  pitch  angle  ^  0°) 
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Figure  A3:  Comparison  of  for  Mach  0.95  and  pitch  angle  =  0° 
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APPENDIX  B 


DETERMINATION  OF  ROLL  BALANCE  BEARING  FRICTION 


VffWWW**  tcTfm 


Determination  of  Rot!  Balance  Bearing  Friction 

The  model  was  removed  from  the  balance  and  replaced  with  a  simple  flywheel.  Compressed 
air  from  an  air  hose  was  then  used  to  spin  the  balance  with  the  flywheel  attached.  Once  the 
balance  had  reached  a  roil  rate  slightly  greater  than  the  maximum  constant  spin  rate  observed 
during  the  wind  tunnel  tests,  the  air  hose  was  removed  and  the  balance  allowed  to  come  to 
rest.  During  this  time,  roll  rate  data  were  recorded  and  saved  to  disk.  A  Hewlett  Packard 
digital  oscilloscope  was  initially  set  to  record  the  expected  high  roll  rates,  and  when  the  roll 
rate  reached  the  lower  end  of  the  timebase  range,  the  timebase  was  changed  so  that  the  lower 
roll  rates  could  be  recorded.  This  procedure  was  carried  out  several  times  and  all  results  were 
combined  into  one  file. 

The  data  for  roll  rate  versus  time  were  differentiated  to  produce  values  for  roll  acceleration 
versus  time.  From  a  knowledge  of  the  inertia  of  the  rotating  parts,  roll  accelerations  were 
then  converted  to  bearing  friction  torques,  to  give  a  set  of  bearing  friction  torques  versus  roll 
rate.  The  graph  of  bearing  friction  torque  versus  roll  rate  is  presented  in  figure  B1  and  has 
been  approximated  in  further  calculations  by  the  two-seginented  straight  line  fit  shown  there. 

This  method  of  measuring  bearing  friction  neglects  some  parameters  (eg.  effects  of  variations 
in  axial  and  radial  loads).  In  this  use,  however,  the  aerodynamic  and  gravity  loads  on  the 
bearings  are  only  a  very  small  fraction  of  the  bearing  design  loads,  and  so  are  assumed  to 
have  a  relatively  small  effect  on  the  bearing  friction  characteristics.  Also,  bearing  friction 
generally  accounts  for  only  a  small  part  of  the  calculated  roll  damping  (in  this  case,  less  than 
5%)  so  that  moderate  uncertainties  in  the  measurement  of  the  bearing  friction  arc  not  a 
significant  contributor  to  the  overall  uncertainty  of  the  roll  damping  measurements. 

In  a  previous  use  of  this  technique,  test  runs  were  conducted  in  an  evacuated  chamber  (SI 
wind  tunnel  working  section)  to  assess  the  contribution  of  aerodynamic  damping.  In  these 
tests  aerodynamic  damping  of  the  flywheel  was  found  to  be  negligible  compared  to  bearing 
friction. 
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APPENDIX  C 
DERIVATION  OF  C 


Derivation  of  Cmq 

The  following  derivation  was  obtained  Mid  modified  from  reference  4. 

Notation: 

a1F  experimental  lift  curve  slope  of  fin 

Cjn^  pitch  damping  coefficient 

Z  reference  length 

d  mortar  diameter 

i  integer  involved  in  summation  to  replace  integral  for  (M,)b 

k,  parameter  in  calculation  of  afterbody  correction  factor  ( 1  -  k») 
kf  correction  factor  for  body  fineness 
1B  body  length 

!□  distance  of  moment  reference  point  aft  of  body  nose 
M  Pitching  moment  about  moment  reference  point 
(M,)b  body  contribution  to  total  pitching  moment  derivative 
(M<,)f  fin  contribution  to  total  pitching  moment  derivative 

q  rate  of  pitch  about  moment  reference  point,  positive  nose  up 
SB  planform  area  of  body 

Sm  cross-sectional  area  of  body  that  contains  maximum  body  width 
Sw  reference  length  (cross-sectional  area  of  mortar  at  maximum  body  width) 

Vc  velocity  of  mortar  relative  to  air  in  undisturbed  flight 
wB  local  body  width 

xB  longitudinal  distance  along  body  axis  to  general  body  station,  measured 
positive  aft  of  moment  reference  point 

xK  longitudinal  distance  along  body  axis  of  V*  chord  point  of  aerodynamic 
mean  clrord  of  fin,  measured  positive  aft  of  moment  reference  point 

p  density  of  air 


In  reference  4  the  pitch  damping  coefficient,  Cny  is  defined  as 

5M/5q 


cm  =■ 


~m<*  Vi  p  Ve  Swc2 

This  equates  to  the  definition  of  pitch  damping,  =  ^j11  ,  which  is  commonly  used  in 

missile  work,  where  two  different  reference  lengtlts  are  used,  namely; 

c  =  d  (missile  diameter)  for  non-dimensionalising  the  pitching  moments 


.  =  d 

and  c  =  ~ 


(missile  radius)  for  non-dimensionalising  the  pitch  damping  coefficient 


> 


r  8M/5q  _  f  6M/Sq  >| 

V4pVeSw(d2/2)  ^pVeSwCl2J 

Cjn^  is  calculated  as  the  sum  of  two  contributions,  namely  body  and  fin,  ie 

Cm,  =  (Hj)b  +  (H))F  . 

Body  contribution,  (M,)b  is  calculated  from 


/2S  \kf(l-k«)  (*‘B’*0  , 

■‘o 

To  perform  the  integral  in  the  equation  above,  the  body  was  divided  into  twenty  transverse 
segments  of  length  1/20;  these  are  numbered  i  =  1  to  20,  from  fore  to  aft.  The  local  body 
width  wBi  and  the  local  moment  arm  xBi  are  determined  at  the  midpoint  of  each  segment.  The 
integral  may  then  be  approximated  as 


ws  Xg  dxB  * 


( 20  \ 
£  wB,  x2 
i  =  1 


Fin  contribution,  (M,)f  is  calculated  from 


(M,)F- -a1F 


d/2 


where  a1F  (in  this  case)  comes  from  wind  tunnel  measured  data. 

Calculated  pitch  damping  coefficients  at  various  Mach  numbers  arc  presented  in  the  table 
below. 


I 


ft 


ft 


ft 


Mach  No. 

(M,)B 

(M,)p 

Cm, 

0.5 

-5.760 

-26.590 

-32.350 

0.7 

-5.926 

-33.808 

-39.734 

0.8 

6.040 

-35.314 

-41.354 

0.9 

6.154 

-43.096 

-49.250 

0.95 

-6.230 

-43.136 

-49.366 

Tabie  Cl :  Pitch  damping  coefficients 


As  seen  in  table  Cl,  the  body  contribution  is  relatively  small  (“15%  of  the  total).  Therefore 
the  uncertainty  in  the  estimated  C™  depends  largely  on  the  quality  of  the  tail  lift  data.  Being 
experimentally  measured,  this  should  be  reliable. 
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